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Phase diagram of NaccCo02 (x ^ 0.71) has been reinvestigated using electrochemically fine tuned 
single crystals. Both phase separation and staging phenomena as a result of sodium mult i- vacancy 
cluster ordering have been found. Phase separation phenomenon is observed in the narrow ranges 
of 0.76 < X < 0.82 and 0.83 < x < 0.86. While x = 0.820 shows A-type antiferromagnetic (A-AF) 
ordering below 22K, x = 0.833 is confirmed to have a magnetic ground state of A-AF ordering below 
~8K and is only reachable through slow cooling. In addition, x = 0.859 is found to be responsible 
for the highest A-AF transition temperature at about 29K. Staging model based on ordered stacking 
of multi- vacancy layers is proposed to explain the hysteretic behavior and A-AF correlation length 
for X - 0.82-0.86. 
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PACS numbers: 74.62. Bf, 74.25. Bt, 74.62. Dh, 74.78. Fk 



I. INTRODUCTION 

Layered material Naa;Co02 has a rich electronic and 
magnetic phase diagram as a function of x, from A- 
type antiferromagnetic ordering for x ^ 0.75, to metal- 
to-insulator transition for x ~ 1/2, and even supercon- 
ductivity for X ^ 1/3 after hydration,^ Although A-type 
AF magnetic ordering transition below 22K has been re- 
ported in all samples of nominal x from 0.75 to 0.85, the 
difference among these concentrations has usually been 
ignored, either due to poorly controlled Na level from 
melt growth or roughly estimated Na content. ^'^^^ The 
high Na vapor loss during high temperature melt growth 
is well known and the diffusive nature of Na ions at 
room temperature makes the control of Na content even 
more difficult,^ which can often lead to an inhomoge- 
neous mixture of phases for x ^ 0.7. Only until recently, 
detailed Na ion ordering has been revealed through neu- 
tron and synchrotron X-ray diffraction studies on single 
crystal samples. The newly found evidence of super- 
structure formed by multi- vacancy clusters in x ~ 0.71 
and 0.84 introduces an idea of doped holes partial local- 
ization, which is able to resolve many intriguing phys- 
ical phenomena found in this layered system, including 
the Curie- Weiss behavior of a metallic system, the en- 
hanced thermoelectric power, the novel spin liquid state, 
the reconstructed Fermi surface, and the origin of A-type 
AF ordering. ''^^^^ Recent studies of x ~ 0.80 and 0.85 by 
Schulze et at. conclude that Na ordering is highly depen- 
dent on the cooling rate, where an additional magnetic 
ordering below 8K appears only after the sample is slowly 
cooled through the 300-200K range. However, the real 
impact of the successive Na rearrangement processes re- 
mains to be clarified and the phase diagram must be 
revisited. 



Herein, using results from additionally improved elec- 
trochemical techniques, specific heat and high resolution 
single crystal synchrotron X-ray Laue diffraction, we re- 
port detailed magnetic and structural phase diagram in 
the region of 0.71 ^ x ^ 0.86. Structural phase separation 
phenomenon is found in the two regions of 0.76-0.82 and 

0. 83-0.86 at room temperature as shown in Fig.[Tl While 
simple hexagonal superstructure of A/TSa is maintained 
in all samples with x in the range of 0.82-0.86, the mag- 
netic ground state turns out to be distinctively different. 
In fact, there are three distinct A-AF transition tem- 
peratures of TAri=22K, TAr2=8K and TAr3=29K found 
in this range, corresponding to a proposed specific multi- 
vacancy layer stack ordering of well-defined stoichiometry 
of X = 0.820, 0.833, and 0.859 respectively, plus x=0.763 
that shows a spin glass like behavior below ~ 3K and with 
a significantly larger superlattice of V^^a. In particular, 
we find x=0.820 to be the most stable phase and cooling 
rate independent, while x=0.833 shows a strong cool- 
ing rate dependent nature with transitions found near 
8K (through slow cooling) and 16K (through fast cool- 
ing). These distinctively different A-AF phases found 
in 0.820, 0.833 and 0.859 can be reasonably constructed 
by adding different levels of additional Na vacancy to 
the ideal di-vcancy formed x=l 1/13=0.846 superstruc- 
ture of V^axvT^ax3cj^i^ Applying a layered staging 
model similar to that used in graphite intercalated com- 
pounds (GlC)fii for example, Tati phase (x=0.820) can 
be described as a stage-2 compound, i.e., where correct 
stoichiometry is obtained by introducing two more Na va- 
cancies into the original ideal V^aXA/T3ax3c super unit 
cell, and these defects create tri-vacancy layers that are 
sandwiched between every two di- vacancy layers. On the 
other hand, TAr2 phase (x=0.833) corresponds to stage-5, 

1. e., tri-vacancy layers are sandwiched between every five 
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FIG. 1: (color online) (a) A revised phase diagram of 
NaxCo02 in the range of 0.71-0.86 and the proposed 
staging models for x=0.820(Tivi=22K), 0.833(TiV2=8K) 
and 0.859 (TiV3=29K). Magnetic ordering temperatures are 
marked by red cross, (b) Na layers with tri- vacancy (T), di- 
vacancy (D) and mono- vacancy (M), where Na2 ions (blue) 
move from the original Na2 site (empty circle) to the Nal 
site (red), (c) Crystal structure for 7-NaccCo02 of PGs/mmc 
symmetry with vacant Nal sites shown in empty circles. 

di- vacancy layers. Most interestingly, the staging model 
suggests that the trivacancy layers serve as nucleation 
centers for the interlayer AF ordering. This picture nat- 
urally explains why x=0.820 has a higher Tn than that 
of x=0.833 because of its shortest inter-trivacancy layer 
distance. The observed phase separation phenomenon is 
a natural consequence of the competing multi-vacancy 
cluster size, super lattice size, and interlayer magnetic 
correlation. 



II. EXPERIMENTAL 

High quality single crystals of well controlled Na con- 
tent were prepared using electrochemical de-intercalation 
technique starting from high Na content crystal of x 
^ 0.84 grown with floating-zone method, the details 
have been documented previously^^i^ii^ The exact Na 
content has been cross checked with c-axis vs. x lin- 
ear relationship constructed from combined high angle 
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FIG. 2: (color online) (a) X-ray diffraction results of (008) 
peak for samples with x in the ranges of (a) 0.83-0.87 and 
(b) 0.76-0.82 at room temperature. The linear x-dependence 
of relative change of (008) diffraction integrated intensities 
suggests the phase separation phenomenon. The broadened 
diffraction peak for x > 0.86 as shown in (a) indicates poor 
ordering. 

TABLE I: Summary of studied NaccCo02 crystal samples 
sample # 1-5 6-10 11-15 16-20 21-25 
X (EPMA) 0.768(4) 0.801(5) 0.820(3) 0.832(2) 0.842(5) 
X (c-axis) 0.771(2) PS 0.820(2) 0.832(2) PS 
X (model) 0.763 PS 0.820 0.833 PS 0.859 

Tiv(K)* 3** 22 22 8 8/29 

*slow cooling **spin glass like behavior 

X-ray diffraction (008) peak position. Inductively Cou- 
pled Plasma (ICP) and Electron Probe Microanalysis 
(EPMA) techniques. ^'^'^ In particular, current study uses 
c(x) linear function that is further calibrated by the phase 
separated boundaries and EPMA is averaged out from 
freshly cleaved crystal surface for more than 100 points. 
A complete list of samples studied is summarized in Ta- 
ble [D Due to the active diffusive nature of Na ions at 
room temperature and the minute differences of Na con- 
tent, all measurements were done on freshly prepared 
crystals within days. Otherwise crystal samples must be 
stored within L-N2 dewar below 200K in order to sup- 
press Na loss from the surface. Synchrotron Laue diffrac- 
tion for Na superstructure investigation was performed 
with synchrotron source in Taiwan NSRRC, and mag- 
netic property characterization was done using Quantum 
Design SQUID MPMS-XL. 



III. RESULTS AND DISCUSSIONS 

While zooming in the region of x > 0.75 using electro- 
chemical technique, we found several concentrations to 
be particularly stable and of two-phase character. As 
indicated in Fig. [21 the evidence of phase separation 
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FIG. 3: (color online) Synchrotron X-ray Laue pattern for 
x=0.801(5) at room temperature, where three twin sets of 
hexagonal units are shown in the reciprocal space (upper in- 
set), which correspond to hexagonal A/TSa (blue), vT2a (red) 
and \/T9a (green) superlattices in the real space (lower inset). 



is demonstrated by the x-dependent evolution of (008) 
diffraction peak integrated intensities, where the growth 
of one end phase is at the expense of the other at the mis- 
cibility gap boundaries without continuous intermediate 
phase in between. As shown in Fig. [1] and Fig. O phase 
separation phenomenon is found to occur in two regions 
of X ~ 0.76-0.82 and 0.83-0.86, where coexistence of two 
end phases grow at the expense of each other as indicated 
by the (008) X-ray diffraction peaks, while solid solu- 
tion connects phases between miscibility gaps. Clearly 
the growth of 0.82 phase is at the expense of 0.76 phase 
for increasing x. Such coexisting pattern of two phases 
can only be observed in the narrow regions of 0.76-0.82 
and 0.83-0.86, although 03-type secondary phase is com- 
monly found for x ^ 0.85 from melt growth.-^*^ Pure 
x~0.86 phase cannot be prepared using electrochemical 
technique, partly due to the fact that it becomes harder 
for Na to intercalate into the host structure of compressed 
c-axisj^i^ Enforcing more Na into the matrix electrochem- 
ically destroys its ordering as indicated by the broadened 
diffraction peaks that correspond to x^O.86 as shown in 
Fig.Ela). 

Fig. [3] shows synchrotron X-ray Laue diffraction re- 
sults for X = 0.801(5), where x sits in the middle of the 
miscibility gap of 0.76-0.82. Although in-house x-ray re- 
sults for X within 0.76-0.82 range all display two-phase 
feature of the two end compounds only, the synchrotron 
X-ray Laue shows there are actually three coexist phases 
at room temperature. As shown in Fig.[3l we can clearly 
identify three twin sets of hexagonal superlattices that 
correspond to A/T9a, vT3a and >/T2a in the real space, 
the latter two can be compared with the published sin- 
gle phase Laue patterns of x = 0.71 (V^a) and 0.84 
(>/T3a).^^^^ The newly found >/T9a superlattice must be 
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FIG. 4: (color online) Magnetic susceptibility measurement 
results for x > 0.80 under applied field of 1 Tesla and 0.01 
Tesla (inset) along the c- and ab-directions respectively. All 
measurements were performed after being zero-field-cooled 
from room temperature through 200K with cooling rate of 
2 K/min. A- type AF ordering is indicated by the cusp of 
c-axis susceptibility Xc and the onsets of low field hysteresis 
shown in the inset. 



due to Na ordering for x ~ 0.76 and has been verified 
by single phase sample, where more Na vacancies are in- 
troduced and a larger superlattice becomes necessary to 
accommodate larger multi- vacancy clusters. In fact it re- 
quires 4.5 vacancies per each superlattice of size A/T9a to 
account for stoichiometry of x ~ 0.76, i.e. x = 1 - ^ 
= 0.763 that is composed of quadri- and pent a- vacancy 
clusters in adjacent layers for 7-NaajCo02^^^ A detailed 
analysis for x ~ 0.76 sample will be presented elsewhere. 
By treating this layered material to be a pseudo-binary 
system composed of Na and C0O2, the triple coexist- 
ing superlattices at constant temperature and pressure 
does not violate the phase rule, in fact it has reached the 
allowed maximum number of three. Considering dom- 
inant domains are from x~0.82 (vT^a with di- vacancy 
clusters) and x~0.76 (A/T9a with quadri/penta- vacancy 
clusters), it's reasonable to have a buffered zone at the 
domain boundary which is built with tri/quadri- vacancy 
clusters of A/T2a superlattice. 

Most of the magnetic susceptibility measurements for 
Naa^Co02 with x ^ 0.75 before show A-type AF signature 
near 22K, i.e., the cusp of Xc under high field. Although 
crude magnetic phase mappings in this range before sug- 
gest that Tat varies between 22-27K,3'i^ Schulze et al 
recently found an additional 8K phase for x ~ 0.80 and 
0.85, which can be obtained only after a slow cooling 
process. With carefully tuned single crystals in the nar- 
row range of 0.82-0.86, we are able to re- visit the mag- 
netic phase diagram and untangle the mystery of Tat 
variation. Magnetic susceptibility measurement results 
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FIG. 5: (color online) Specific heat measurement results for x 
= 0.820(3) and 0.832(2) obtained from slowest cooling rate of 
0.5 K/min. The inset shows x=0.832(2) has a strong cooling 
rate dependence, where 8K phase grows at the expense of 16K 
phase as cooling rate reduces from 15 to 0.5 K/min. 



are shown in Fig. [H all measurements were done after a 
slow cooling rate of 2 K/min. We find that Tn does not 
change with x monotonously and continuously, instead, 
the four phases at the two miscibility gap boundaries 
shown in Figs. [1] and [4] are responsible for the different 
characteristic Tat's, where x ~ 0.76 shows Curie- Weiss 
behavior down to 5K (not shown). The onsets of A-AF 
transitions are indicated by the cusp of Xc under high 
filed, which occur at TAri=22K, TAr2=8K and TAr3=29K 
for X = 0.820, 0.833 and 0.859 respectively, while 0.801(5) 
and 0.842(5) data reflect their mixed phase nature, i.e., 
superposition of contributions from the end members 
of 0.76-0.82 and 0.83-0.86 miscibility gaps respectively. 
There is ZFC/FC irreversibility found below Tn for both 
Xc and Xab at low field, although stronger FM satura- 
tion moments are seen along the ab-direction. Such A- 
AF ordering has been verified by the neutron scatter- 
ing for Nao.82Co02,^ where strong field dependence of 
magnetization below Tn ^ 22K has been confirmed to 
be metamagnetic.^^ Current low field measurement is in 
agreement with that reported for x ~ 0.85,^^ although 
our data indicate that the 8K phase is coming from the 
phase of x closer to 0.833, while a more stable phase of 
22K is from x closer to 0.820. 

We find that the different onsets of A-AF transition 
between x = 0.820(3) and 0.832(2) are clearly demon- 
strated by the cooling rate dependence of Tat as revealed 
by the specific heat data shown in Fig.[5l where 22K tran- 
sition for x=0. 820(3) is independent of cooling rate, while 
fast cooling rate moves Tn discretely from 8K to 16K for 
x=0. 832(2). 16K phase occurs in x=0. 832(2) when a fast 
cooling is applied, while it decreases at the expense of 



8K phase generation under decreasing cooling rate, al- 
though the existence of minor 22K is difficult to avoid 
completely for samples near the phase separated bound- 
ary. The ratio of the minor 22 K phase to the major 8 K 
phase in the present x=0.833 sample can be estimated to 
be 3.6% from the entropy change of the small anomaly at 
22 K for Nao.833Co02. Sodium ion diffusion is active at 
room temperature for high x samples,^ but it freezes be- 
low ~200K as indicated by the sharp increase of 1/Ti for 
^^Na due to Na motion.^^ Sample of x=0. 832(2) must be 
cooled through the temperature range of 300-200K with a 
rate slower than 10 K/min in order to reach the magnetic 
ground state that corresponds to 8K magnetic ordering. 

The entropy associated with the 22 K transition in 
Nao.82oCo02 and the 8 K transition in Nao.833Co02 
are estimated to be AS~170 mJ/mol K and AS~200 
mJ/mol K, respectively. Taking the surface value of esti- 
mated AS and the transition width character, these re- 
sults indicate that the magnetic moments order better in 
Nao.833Co02 than in Nao.82oCo02. The relatively poor 
ordering of stage-2 for x=0.820 (TAri=22K) than that 
of stage-5 x=0.833 (TAr2=8K) is interestingly in agree- 
ment with the alternating T-Q stacking requirement as 
described by the x=0.71 superstructure model before^ii^ 
i.e., tri- vacancy is not favorable in the even-layer within 
P63/mmc symmetry and there must exist mixing stages 
of 1 and 3 for x=0.820. The AS value for x=0.833 is 
more than 10 times larger than that reported in Ref. [10], 
which suggests a nearly single 8 K phase in the present 
sample. On the other hand, these measured AS values 
are only about 20% of the entropy estimated from the 
complete ordering of fully localized spin-1/2 Co^+ ions. 
This discrepancy might indicate the failure of the simple 
ionic Co^+-Co^+ picture. 

Since all samples with x ^ 0.82 show identical su- 
perlattice size of V^^a, the subtle difference for x ~ 
0.82 and 0.83 must be related to the Na rearrangement 
generated by the additionally introduced Na defect that 
causes deviation from the ideal di-vacancy constructed 
0.846=11/13 of V^axyT3ax3c superstructure.^^^^ But 
what kind of mechanism is responsible for these discrete 
Tat's of Ax only 1-3% apart? The secret lies in the stack 
ordering of 2D hexagonal superlattices. From our previ- 
ous studies on the structure of 0.71 and 0.84,^^^^ the ideal 
superlattice has a 3c periodicity. The 3c periodicity for 
x=0.820 is once again confirmed by electron diffraction 
patterns on single domain crystals as shown in Fig. [6l 
Although [001]p diffraction pattern cannot tell the peri- 
odicity along c-direction, perfect indexing for diffraction 
patterns with transmitted beam along primitive [101]p 
and [201]p can only be achieved with the help of 3c peri- 
odicity assignment. When one and two more Na defects 
per 3c unit (i.e., six layers of Na) are introduced into the 
perfectly ordered original 0.846 = 11/13 superstructure,— 
two additional stoichiometrics of 0.833 = 0.846— | x ^ 
and 0.820 = 0.846— | x ^ can be introduced, as verified 
by our X-ray and magnetic measurement results shown 
above. When magnetic ordering occurs, spins from itiner- 
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FIG. 6: (color online) ) Electron diffraction patterns of sin- 
gle domain Nao.82oCo02 for transmitted e-beam along [001]p, 
[101]p, and [201]p projections (p, denoting the primitive lat- 
tice). The superlattice spots surrounding the intense primi- 
tive cell reflections in the ab-plane projection, [001]p, indicate 
the superlattice ordering of vT^aXA/TSa without c-axis infor- 
mation. In the azimuthal projections of [101]p and [201]p, the 
superlattice reflections suggests 3c ordering, i.e. can only be 
indexed correctly with V^^ax V^^axSc superlattice. 



origin of the ferromagnetic layers which then order anti- 
ferromagnetically between layers to form the A type AF 
ordering. Here we suggest that the driving force for in- 
ter layer coupling may lie in the T-layer. The tri- vacancy 
has one extra negative charge which lowers the tunneling 
barrier between the localized holes on the adjacent layers 
and enhances the antiferromagnetic spin correlation be- 
tween them. Thus the T-layer may form the nucleation 
layers to drive the three dimensional AF order. This pic- 
ture explains why Tati is 22K for x=0.820 vs. is 
8K for x=0.833 where the spacing between T-layers is 
much larger. Upon rapid cooling, some stage 4 and stage 
6 states may form. The stage 4 meta-stable phase may 
be responsible for the intermediate Tat of 16K. The hys- 
teretic behavior observed below Tn (see Fig. |4]) could 
also be explained by either the in-plane FM domain ef- 
fect or by uncanceled A-type AF moments along the c- 
axis as a result of mixed staging. The phase separation 
observed near 0.83-0.86 can also be explained using the 
same stage model. Since 0.859=0.846+^ x i.e., one 
more Na ion (not vacancy) is introduced into the original 
ideal x=0.846 phase of ^/TSax A/T3ax3c superstructure, 
the di- vacancy is converted to a mono- vacancy (M) form- 
ing a stage 5 stacking. The hole density on either side of 
the M-layer is now reduced by ^ x ^ and we may expect 
an even strong tendency toward Stoner ferromagnetism. 
The higher transition temperature of the Co layers adja- 
cent to the M-layer may explain the higher Tats = 29K 
for x=0.859. 



ant electron or localized electrons near Co ions are certain 
to be affected by the rearrangement of Na multi- vacancy 
clusters in the nearby layers. In fact the in-plane inter- 
vacancy cluster distance ^/TSa is nearly twice the inter 
C0O2 distance. Since every one extra Na vacancy in- 
troduced into the ideal 3c unit would convert the orig- 
inal layer of di-vacancy formed 2D superstructure into 
tri-vacancy, we can thus simplify the stacking problem 
into stack ordering between the di-vacancy (D) and tri- 
vacancy (T) layers along the c-direction. We propose a 
new staging model as shown in Fig. [1] to explain such 
stack ordering, which shows strong resemblance to the 
staging phenomenon often observed in the 2D interca- 
lated graphite compounds. The x=0.833 phase which 
has only one Na defect introduced could have a stage-5 
construction, while x=0.820 phase of two Na defects per 
3c unit must have a stage-2 construction, i.e., there are 
five and two D-layers sandwiched in between T-layers. 

We can now use the staging picture to interpret the 
variety of magnetic ordering temperatures observed in 
the range 0.82 ^ x ^ 0.86. As discussed previously,'' 
the di-vacancy may localize a carrier on the adjacent Co 
layer, leaving a low density hole gas of density ^, which 
is unstable to Stoner ferromagnetism. This may be the 



IV. CONCLUSIONS 

In conclusion, we have revised NaajCo02 phase dia- 
gram in the range of 0.71-0.86 using electrochemically 
fine tuned single crystal samples. The puzzling and in- 
consistent measurement results in this range before have 
been clarified and interpreted as a result of phase separa- 
tion and staging phenomena. A direct link between the 
high temperature Na ion (vacancy) ordering and the low 
temperature magnetic properties have been established. 
The highly correlated ion, magnetic and charge orderings 
in layered Naa;Co02 can provide invaluable information 
to the study of strongly correlated electron low dimen- 
sional system which has itinerant electrons on a triangu- 
lar lattice. 
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